Precessing vortices and antivortices in ferromagnetic elements 
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A micromagnetic numerical study of the precessional motion of the vortex and antivortex states 
in soft ferromagnetic circular nanodots is presented using Landau-Lifshitz-Gilbert dynamics. For 
sufficiently small dot thickness and diameter, the vortex state is metastable and spirals toward the 
center of the dot when its initial displacement is smaller than a critical value. Otherwise, the vor- 
tex spirals away from the center and eventually exits the dot which remains in a state of in-plane 
magnetization (ground state). In contrast, the antivortex is always unstable and performs damped 
precession resulting in annihilation at the dot circumference. The vortex and antivortex frequencies 
of precession are compared with the response expected on the basis of Thiele's theory of collective 
coordinates. We also calculate the vortex restoring force with an explicit account of the magne- 
tostatic and exchange interaction on the basis of the 'rigid' vortex and 'two-vortices side charges 
free' models and show that neither model explains the vortex translation mode eigenfrequency for 
nanodots of sufficiently small size. 

PACS numbers: 75.70.Kw,75.75.Fk,75.75.Jn,75.78.Cd,75.78.Fg 
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I. INTRODUCTION 

The vortex state is one of the equilibrium states of thin 
soft ferromagnetic elements of micrometer size and below 
(magnetic dots). The interplay between the magneto- 
static and exchange energy favours an in-plane, closed 
flux domain structure with a 10 — 20 nm central core, 
where the magnetization turns out of plane to avoid the 
high energetic cost of anti-aligned moments. Core rever- 
sal can be triggered by application of an in-plane pulsed 
field or pulsed current allowing the possibility of applica- 
tion of patterned thin film elements in data storage and 
magnetic and magneto-electronic random access mem- 
ory [Ij. Core reversal is usually assumed to arise from 
the spontaneous creation of a vortex-antivortex (VA) pair 
(vortex dipole) of opposite polarity with respect to the 
original vortex, followed by collision of the pair with the 
original vortex. A fundamental understanding of the dy- 
namics of vortices and antivortices is therefore necessary 
to control the switching of the magnetization. 

The basic excitation mode of the vortex or antivortex 
state from its equilibrium position is in-plane gyrotropic 
motion. It is a low frequency (GHz) mode corresponding 
to the displacement of the whole structure. The gener- 
alized dynamic force can be determined using Thiele's 
collective- variable approach [2]-[3]- Theoretical [1]-[S] 
and experimental [6] studies of the dynamics of magnetic 
vortices in 2D films have shown a connection with the 
topology of the magnetization structure . Magnetic vor- 
tices confined in circular dots can be described by ana- 
lytical models based on different methods for accounting 
for the magnetostatic interaction [7]-[S]. The vortex may 
be 'rigid' or deform so that no magnetic charges appear 
at the side of the cylinder. The latter (two- vortices side 
charge free model) provides a good description of the dy- 
namic behavior of vortices in submicron-sized permalloy 



dots, in particular the increase of vortex eigenfrequency 
with dot aspect ratio L/R, where L is the dot thickness 
and R is the dot radius. The basic assumption in these 
calculations is that the vortex displacement I from equi- 
librium, at the dot center, is small I << R. 

The main objective here is to consider the response of 
the vortex state to large perturbations, induced for in- 
stance by thermal activation. We focus our attention to 
dots of sufficiently small radius and thickness so that the 
vortex state is metastable . Using micromagnetic calcu- 
lations it is shown that vortex stability can be defined 
in terms of a critical displacement Ic leading to an irre- 
versible transition to the ground state characterized by 
in-plane magnetization. We test the accuracy of the col- 
lective coordinate representation for vortex and antivor- 
tex dynamics and discuss the limitations of the approxi- 
mate analytical models [8],[9j. 



II. VORTEX AND ANTIVORTEX PRECESSION 

The dynamical behavior of a single vortex or antivor- 
tex trapped in a ferromagnetic nanodot was studied 
using a finite-difference Landau-Lifshitz-Gilbert micro- 
magnetic model. The material parameters for permal- 
loy (Ni8oFe2o) were used in the calculations: saturation 
magnetization Alg — 800 emu/cc and exchange stiff- 
ness coefficient A = 1.3 x 10~^ erg/cm. The mag- 
netic anis otropy was neglected and the exchange length 
is lex = \/A/(27rM2) = 5.7 nm. The disk thickness is in 
the range L = 5 — 10 nm , comparable to the exchange 
length, so the magnetization dependence along the dot 
normal can be neglected to a first approximation. The 
disk radius is not large compared to the vortex core and 
is initially taken as i? = 30 nm. Inspection of the phase 
diagram of magnetic ground states for this material [10] 
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FIG. 1: a) Vortex snd b) antivortex structure in permalloy 
dot of thickness L = 5 nm and radius _R = 30 nm. 



indicates that for the particular choice of disk dimensions, 
in-plane magnetization is the ground state of the system. 
The integration of the LLG equation is carried out using 
a damping parameter a = 0.01, which is appropriate for 
permalloy pT| . 

The initial vortex or antivortex structure is defined as 
follows: Using the coordinate system with the z axis par- 
allel to the dot cylindrical axis, the magnetization compo- 
nents are tti^ = \cosQ{p)^mx + irriy = sinOexp(iK((/) — 
0o)), where m(r) = M(r)/Afs is the reduced magneti- 
zation, p, (f) are cylindrical coordinates, A = ±1 is the 
vortex polarity, and Q{p) is the magnetization angle to 
the normal (z) direction. The vortex number is k = 1 
for a vortex and k = — 1 for an antivortex structure. 
The constant (j>o = ±7r/2 defines the chirality i.e. the 
direction of the curl of the vortex. The specific choice 
of chirality is largely insignificant for vortex dynamics, 
in contrast to the vortex number k and polarity A which 
play an important role. 

The initial profile 8(p) is assumed to have the model 



form COS0 — {cosh{cp/lex))~^ which is appropriate for 
easy-plane ferromagnets ^12, . The parameter c defines 
the core radius b. The core radius is determined by the 
dot thickness L and the exchange length l^x but is inde- 
pendent of the dot radius 13^ . In the following, dot thick- 
ness is i = 5 nm, unless specified otherwise. Numerical 
calculations of relaxed vortex and antivortex structures 
at the dot center were found to be in good agreement 
with the model profile for c ~ 1. Typical examples of a 
vortex and antivortex structure are shown in Fig. [l] It 
should be noted that other choices for the vortex pro- 
file are possible [2], for instance tan(9/2) = p/b was 
obtained by Usov et al. |15j using a variational proce- 
dure to minimize the exchange energy whereas the vor- 
tex core radius b was determined by minimization also 
of the magnetostatic energy. This form is employed in 
analytical calculations of vortex states in ferromagnetic 
disks [2]-[9], however, micromagnetic modelling by Scholz 
et al. [lOj . confirmed by our calculations, has shown that 
it underestimates the core radius b, defined in Ref. [lU] as 
the first rjiz = crossover from the vortex center. 

The dynamical behavior of the magnetic vortex (k = 1) 
depends on the initial displacement l{t = 0) = lo from 
the disk center. If the displacement is equal or smaller 
than some critical value Ic = 0.52R, the damped preces- 
sion of the vortex leads to relaxation to the disk center. 
For positive vortex polarity A = 1, an anticlockwise pre- 
cession is observed as shown in Fig. [2^. At the start of 
the simulation the precession is not smooth as a result of 
the internal relaxation (deformation) of the model vortex 
to minimize the energy during precession. If lo > Ic, the 
damped precession is clockwise and the distance from the 
dot center increases, as shown in Fig. [2]d, until the vor- 
tex is annihilated and the magnetization is aligned along 
the in-plane direction with quasi-uniform magnetization, 
the so-called 'leaf state [TB] which is the ground state 
of the system. Irregularities in the precessional motion 
arise from the uncertainty on the position of the vortex. 

The antivortex instead is always unstable. For any 
choice of initial displacement ^o, the antivortex performs 
damped precession to the edge of the disk and is an- 
nililated. For positive polarity, anticlockwise precession 
is observed. The sense of gyrotropic motion of a vortex 
or antivortex is switched on reversal of the polarity. 



III. THE COLLECTIVE COORDINATE 
APPROACH 

The damped precession of the vortex or antivortex can 
be described using Thiele's equation [2],|5 augmented by 
a dissipative term. 
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where 1 = {lx,ly) is the position of the vortex center 
and E{lx,ly) is the potential energy of the shifted vor- 
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I = ^/l"^ + ly , dE/dlx = E'lx/l and it is straightforward 



(2) 



to show that 

Ix — 2Qrily = -uoly 

2Qr]lx + ly = iolx 

where the angular frequency is 

1 BE 
'^~ Gl~dl 

The vortex motion in complex form is given by 

(1 + 2Q'qi) \lx+ ilyj = iijJ {Ix + Hy) 
Introducing polar coordinates Ix + ily = le^'^ 
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FIG. 2: Trajectory of a vortex of positive polarity A = 1 in 
zero field for a time interval < i < 2 x 10''. The initial 
position vector of the vortex is a) lo — (0.527?, 0) and b) 
(0.53/?, 0). The anisotropy is neglected and the damping is 
a = 0.01. 



tex. The gyroforce G x dl/dt depends on the topolog- 
ical structure of the magnetization and is proportional 
to the gyrovector G = — Gz, where the gyroconstant is 
G = 2TTn\LMsh and 7 = 1.76 x lO'^ rad Oe~i s~^ is the 
gyromagnetic ratio. Q = —huX is the skyrmion number 
and 77 is the dissipation constant. 

For axially symmetric energy potential E = E{1) where 



^(0 

1+7^2 



(7) 



Dividing and integrating over the time interval of the 
damped precession, the time dependence of the preces- 
sion angle is 



0(i) 



V Vo 



(8) 



where lo is the initial displacement of the vortex center. 
The clockwise or anticlockwise sense of gyration is there- 
fore dependent on the skyrmion number Q. 

Micromagnetic simulations of vortex motion were per- 
formed and the position of the vortex center (Ix, ly) was 
determined by a method of interpolation for the position 
of maximum ttIj,. The precession angle = arctan(/j,//2;) 
was found to vary linearly with the logarithm of the vor- 
tex shift, in agreement with Eq.(l8]). Fig. [s] shows numeri- 
cal data for damped precession of a vortex of positive po- 
larity {Q = —1/2) with initial displacement lo — 0.52R. 
The relaxation to the disk center involves many revolu- 
tions (Fig. [2^) and the gradient — 1/ry provides an accu- 
rate estimate of the dissipation constant t] = 0.013. It is 
evident that for permalloy nanodots, the damping in the 
vortex motion is weak and the angular frequency of pre- 
cession in Eq.(l7]) can be approximated using uj ~ d(f>/dt. 

Fitting the numerical 4>(t) curve to a 4th order poly- 
nomial, the time variation of the angular frequency uj{t) 
can be determined. The vortex shift l{t) exhibits oscilla- 
tions that are neglected by fitting to a 4th order polyno- 
mial and the uj{l) dependence is obtained using the uj{t) 
curve. The same procedure is employed for damped pre- 
cession leading to vortex annihilation {lo = 0.53i?). The 
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FIG. 3: Azimuthal angle of the vortex position as a function 
of the logarithm of the reduced vortex displacement s = l/R. 
The initial radial position is lo = 0.52R. 
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FIG. 4: Precession frequency of a vortex (jj/2n as a function 
of the reduced diplacement s — l/R from the center of the 
dot. 



combined results for the dependence of the frequency of 
precession w/27r on reduced vortex shift s — l/R arc il- 
lustrated in Fig. |4] 

For small displacement of the vortex center from its 
equilibrium position (l — 0), the potential energy is 
E{1) = E{0) + {\/2)kP, where k is the stiffness coeffi- 
cient and the eigenfrequency is Uo = n/G 0. At the 
critical displacement Ic, corresponding to a maximum in 
the potential energy E{1) the precession frequency van- 
ishes. 
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FIG. 5: (a) Reduced energy density of a permalloy dot 
e = E /A-kM'^V as a function of normalized vortex displace- 
ment s. Results are shown for two sets of micromagnetic simu- 
lations (markers) and analytical (solid lin e) calculations using 
Thiele's collective variable theory (Eq |10[ ). The maximum en- 
ergy density emax occurs at vortex displacement Sc- (b) The 
contribution of the magnetostatic and exchange terms to the 
total energy, obtained from micromagnetic calculations. 



The motion of vortices and antivortices is driven by 
the restoring force dE/dl (Eq[T]). The potential energy 
of the shifted vortex is axially symmetric E = E{1) and 
can be written 



E{1) = ^(0) + G / uj{p)pdp 



(9) 



A simpler form in terms of the reduced energy density 
e = E/At:M'^V over the dot volume V is 



e{l) = e(0) - 4g / uj{s')s'ds' 
Jo 



(10) 



where s' = p/R, s = l/R is the reduced displacement 
of the vortex and the time associated with uj is scaled 
by To = l/AiT-fMs. The function uj{s') was determined 
from micromagnetic simulations of the vortex precession 
and using Eq([lO| the curve e{l) expected from appli- 
cation of the theory of collective coordinates, was ob- 
tained, as shown in Fig. [5k. Superimposed is the energy 
evaluated directly from simulations of vortex relaxation. 
Thiele's theory appears to provide a good description of 
the vortex precession despite its limitations, for instance, 
a) Thiele's approach is known to be an approximate de- 
scription of vortex dynamics in infinite films b) the vortex 
is here confined in a nanodot {R = 5.Mex) and c) the vor- 
tex structure does not remain rigid during the relaxation 
process but is modified as a result of the change in the 
distribution of the demagnetizing fields. The magneto- 
static and exchange contribution to the energy variation 
e{l) was obtained from micromagnetic calculations and 
is shown in Fig. [5]3. It should be noted that incorporat- 
ing the demagnetizing energy to the total anisotropy, is 
strictly valid for infinite thin films and results in mono- 
tonically decreasing energy e(/), as reported in Ref. [T7] . 
The oscillations in the energy variation during precession 
are related to the finite micromagnetic grid. 



The potential energy of the vortex attains a maximum 
value at some critical value of the displacement Sc = 0.52 
corresponding to a zero crossover of the precession fre- 
quency w (Fig. |4| . The stability of the vortex at the dot 
center arises from the magnetostatic energy, in particular 
the volume magnetic charges resulting from vortex defor- 
mation and the surface charges at the side of the cylinder 
[9] . The face charges do not depend on s since the charge 
distribution on the top and bottom surfaces of the disk is 
unchanged with the vortex displacement. The exchange 
energy decreases with increasing vortex shift s [S]. The 
magnetostatic and exchange contributions to the restor- 
ing force are in exact balance at the point of maximum 
energy. 



A similar analysis was carried out for an antivortex 
structure in a dot of identical dimensions. The potential 
energy decreases monotonically with increasing displace- 
ment s, as shown in Fig. [6] Application of the collec- 
tive coordinates treatment results in the solid curve in 
Fig. |6] of slightly smaller curvature. The vortex insta- 
bility arises from the uncompensated magnetic charge 
distribution within the antivortex core (Fig. IIId) , so the 
magnetostatic energy is reduced upon motion away from 
the dot center. The precession frequency increases dur- 
ing the relaxation process as shown in Fig. [7] as a result 
of the steeper energy gradient for large displacement s. 
Assuming identical position, it is evident that the an- 
tivortex precesses faster than the vortex as a result of 
the larger magnetostatic energy gradient. 
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FIG. 6; Energy density e — E/iivM^V vs antivortex displace- 
ment s. The notation is similar to Fig. [Sk. 
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FIG. 7; Antivortex precession frequency (jj/2n as a function of 
the reduced diplacement s — l/R from unstable equilibrium 
position at the dot center. 



IV. DEPENDENCE OF VORTEX PRECESSION 
ON DISK SIZE 

The maximum potential energy E^ax of the shifted 
magnetic vortex, evaluated from plots such as Fig. [5k, 
depends on the radius of the dot, where it is confined. 
In Fig. [8] micromagnetic calculations of the reduced dot 
energy density e are shown as function of dot radius R, 
scaled by the exchange length. The curves correspond 
to the maximum vortex energy Emax and the minima 
associated with the three equilibrium states of the mag- 
netization (in-plane, perpendicular, vortex). The vortex 
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FIG. 8: Micromagnetic calculations of scaled dot energy den- 
sity e vs dot radius R (in units of lex) for the three equilibrium 
states of the magnetization and the vortex state of maximum 
energy. 
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Barrier to vortex escape (es), defined by the relation 
jortex as a function of dot radius (in units of 



FIG. 10: Reduced displacement for maximum energy of the 
vortex state Sc as a function of dot radius (in units of lex)- 



critical size for vortex instability is larger than Rg since 
the latter is defined assuming a random perturbation dif- 
ferent than shifting the whole vortex. The vortex barrier 
energy increases for larger dots and attains a maximum 
value, for R > lOlex, related to the vortex annihilation 
field [7]. The displacement Ic for vortex escape increases 
with disk size to the maximum value imposed by the 
disk perimeter {Ic/R — >■ 1), attained for sub-micron dots 
R» lex- For sufficiently large dots, the vortex is within 
the domain of attraction of the dot center, irrespective 
of the initial position. 

Micromagnetic calculations of the dependence of the 
fundamental vortex eigenfrequency uJo/^tt, obtained for 
small perturbation s « 1 as in F ig.jil on the dot aspect 



ratio P = LjR are shown in Fig. 11 The dot thickness 



was fixed {L = 10 nm) and the radius R was allowed to 
vary. The eigenfrequency attains a maximum value at 
/3 = 0.35 and vanishes for smaller radius /3 ~ 0.5 when 
the vortex becomes unstable. The maximum value arises 
from the change in the relative contribution of the magne- 
tostatic and exchange energy to the stiffness coefficient 
/i. For instance, the eigenfrequency assuming a 'rigid' 
vortex is [S] 



state is unstable for small dots R < R^, metastable for 
Rs < R < Req and a ground state for R > R^q where 



the values Rs 



2.51 



ex and Req 



16lex are obtained 



for the absolute and equilibrium single domain radius re- 
spectively. The corresponding variation of vortex barrier 
energy es = Emax - E^ortex and displacement U = ScR 
for vortex escape are shown in Figs.[9]and 10 respectively. 
For small dots R = i-blex , the vortex is unstable and any 
shift from equilibrium at the dot center results in relax- 
ation to the ground state (in-plane magnetization). The 
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where Fi{x) — j dtt^^ f{xt)Ji{t) corresponds to the av- 
eraged in-plane dot demagnetizing factor, f{x) = 1 — 
[1 — exp(— a;)]/a; and Ji is the Bessel function. Using 
the approximation i^i(/3) ~ (/3/27r)[ln(8//3)] - 1/2], valid 
for /3 << 1, it can be shown that the eigenfrequency is 
maximum at radius R = Air / L[\n{8 / /3) — 3/2]~^ where 
all lengths are in units of lex- Previous studies were re- 
stricted to sub-micron sized dots where the second term 
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FIG. 11: Micromagnetic (markers) and analytical (solid lines) 
calculations of vortex eigenfrequency Wo/27r as a function of 
dot aspect ratio /3 — L/R a) the rigid vortex model and b) 
the two-vortices model. The dot thickness is here L = 10 nm. 



of Eq.( 11 ) could be neglected [9], so a monotonically in- 
creasing eigenfrequency Wo(/3) arising from the magneto- 
static energy only was reported. 

The micromagnetic calculations are compared in 
Fig. [Tl] with the curves obtained using the 'rigid' vor- 
tex and 'two-vortices' approximations for the magneto- 
static energy. For sub-micron sized disks with aspect 
ratio /3 < 0.05 corresponding to a radius R > 200 nm, 
the micromagnetic calculations are in good quantitative 
agreement with the 'two vortices' model, as reported in 
Ref. [9 . In this regime, the eigenfrequency is determined 
primarily by the magnetostatic energy. The rigid vor- 
tex approximation fails to describe the dynamic behav- 
ior since the magnetostatic energy can be decreased by 
elimination of the surface charges at the disk perime- 
ter at the expense of some contribution from volume 
magnetic charges arising from vortex deformation. For 



smaller disks R < 200 nm, the vortex eigenfrequency is 
between the predictions of the two models. A reduction 
in side charges occurs but is not complete as a result of 
the large expense in exchange energy arising from vortex 
deformation. Similar results can be obtained in principle 
for thinner disks, however, the 'two vortices' approxima- 
tion is then valid for larger cylinders which are not easily 
amenable to micromagnetic simulations. 



V. CONCLUSIONS 

Micromagnetic calculations were carried out of the pre- 
cessional behaviour of a single magnetic vortex or an- 
tivortex confined in a permalloy circular nanodot. The 
existence of two domains of attraction for the vortex state 
are identified arising from a maximum in the potential 
energy of the shifted vortex. This effect is atrributed 
to the competition between the magnetostatic attractive 
and exchange repulsive forces on the shifted vortex. An- 
tivortices instead are always unstable and trace a spi- 
ral trajectory of increasing distance from the dot center 
followed by annihilation at the dot envelope. The pre- 
cessional behaviour of vortices and antivortices is satis- 
factorily described by Thiele's theory of collective coor- 
dinates, relating the angular frequency of precession to 
the gradient of the potential energy. For small nano- 
sized dots, however, the 'rigid' vortex and 'two-vortices' 
approximation for the magnetostatic energy is not sat- 
isfactory and the exchange forces have a significant ef- 
fect on the translation mode vortex eigenfrequency. For 
antivortices, the development of an analytical model to 
account for the magnetostatic interaction in circular dots 
is clearly needed to provide further insight on the results 
of our micromagnetic calculations. 

Vortex stability is necessary in applications of nanos- 
tructured patterned media for data storage. Microfabri- 
cation downscaling implies that the displacement Ic and 
associated energy barrier may become useful characteri- 
zation parameters of the thermal stability of the recorded 
information. 
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